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Abstract 


Interdisciplinary concept of dissipative soliton is unfolded in connection with ultrafast fi- 
bre lasers. The different mode-locking techniques as well as experimental realizations of 
dissipative soliton fibre lasers are surveyed briefly with an emphasis on their energy scal- 
ability. Basic topics of the dissipative soliton theory are elucidated in connection with 
concepts of energy scalability and stability. It is shown that the parametric space of dissi- 
pative soliton has reduced dimension and comparatively simple structure that simplifies 
the analysis and optimization of ultrafast fibre lasers. The main destabilization scenarios 
are described and the limits of energy scalability are connected with impact of optical tur- 
bulence and stimulated Raman scattering. The fast and slow dynamics of vector dissipa- 
tive solitons are exposed. 


Keywords: Ultrafast fibre laser, mode-locking, dissipative soliton, non-linear dynamics, 
vector solitons, optical turbulence, stimulated Raman scattering 


1. Introduction 


Over the last decades, ultrafast fibre laser technologies have demonstrated a remarkable 
progress. By definition [1-4], these technologies concern generation, manipulation and 
application of optical pulses from a fibre laser or a laser-amplifier system with (i) peak power 
P, exceeding substantially an average laser power Pw and (ii) pulse widths t which are much 
lesser than a laser round-trip period T. Such a definition can be re-interpreted in terms of a 
laser mode-locking, which means a phase-locked interference of the M laser eigenmodes 
producing an equidistant train of ultrashort pulses. Then, t œ 1/Môw (M is a number of locked 
eigenmodes and ôw is an inter-mode frequency interval defined by T œ 1/6w) and Py œ MP me 
[5,6]. It means that a mode-locked laser generates a comb of equidistant optical frequencies 
comprising the broad spectral range A oc Mow. It is clear that the substantial enhancement of 
P, (by the factor of M œ 1/tôw ~10° + 10°, i.e. upto over-MW level [7,8]) and t-reduction (~T/M, 
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i.e. down to sub-100 fs level [9]) promise an outlook for different applications [10] including 
non-linear and ultrasensitive laser spectroscopy [11-14], biomedical applications [15-19], 
micromachining [20,21], high-speed communication systems [22], metrology [23] and many 
others [24,25]. The extraordinary peak powers in combination with the drastic pulse width 
decrease bring the high-field physics on tabletops of a mid-level university lab [26-29]. 
Moreover, the over-MHz pulse repetition rates dw provide the signal rate improvement factor 
of 10° + 10* in comparison with that of classical chirped-pulse amplifiers [26]. As a result, the 
signal-to-noise ratio enhances substantially, as well. 


Another aspect of the ultrafast laser applications is connected with studying non-linear 
phenomena [30]. Ultrafast lasers became an effective platform for investigation of general non- 
linear processes such as instabilities and rogue waves [31,32], self-similarity [33] and turbu- 
lence [34]. A coherent self-organization in such non-linear systems [35,36] is the keystone of 
this review, and it will be considered below in detail. But here, we have to point at the 
multidisciplinary context of our topic. The ultrafast fibre lasers can be treated as an ideal 
playground for exploring of non-linear system phenomenology as a whole [37]. Such a 
playground spans gravity and cosmology [38], condensed-matter physics and quantum field 
theory [39-41], biology, neurosciences and informatics [42,43]. The advance of ultrafast laser 
technology is that the theoretical insights promise to become directly testable, controllable and, 
on the other part, the theory can be urged by new precise measurable experimental challenges. 


To date, the solid-state lasers allowed generating shortest pulses with highest peak powers 
directly from an oscillator with high repetition rates (6w > 1 MHz) [44-49]. The main advan- 
tages of solid-state laser systems are (i) broad gain bands (i.e. very large M) allowing generation 
of extremely short pulses (t approaches one optical cycle for such media as Ti:Sp or Cr:chal- 
cogenides) and (ii) covering the spectral range from visible (Ti:Sp) through infra-red (Ti:Sp, 
Cr:forsterite and Cr:YAG) to mid-infrared (Cr:chalcogenides) wavelengths, as well as (iii) 
possibility of independent and precise dispersion [50] and non-linearity [46,48] control. 
Nevertheless, fibre lasers have unprecedented prospects [51] due to (i) possibility of mean 
power scaling provided by large gain, (ii) high quality of laser mode, (iii) reduced thermo- and 
environment-sensitivity, (iv) compactness and integrity of laser setup. Additionally, one has 
to point at broader gain bands of fibre media in comparison with the energy-scalable, thin- 
disk, solid-state oscillators operating within analogous wavelength ranges [9,52]' and possi- 
bility to break into deep-UV and mid-IR optical spectral ranges [55,56]. 


In this review, we will concern the concepts of a mode-locking and a dissipative soliton in a 
nutshell. 


2. Mode-Locking 


The concept of mode-locking is universal and closely connected with a principle of synchro- 
nisation of coupled oscillators [57—64]. A laser is, in fact, the interferometer which possesses a 


1 The width of a gain band is not a decisive factor per se because both pulse width and its spectrum are affected by various 
factors including higher-order dispersions, non-linearity, etc. [53,54]. 
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set of eigenmodes (longitudinal modes) separated by ôw = 27/T. Simultaneously, it is an active 
resonator, which means an amplification AA of mode amplitude A during the resonator round- 
trip in the vicinity of the maximum gain frequency ay, as 
AA(w)=(g (w) -?) A(@) -a(w -w)2A(w), where g(w,) is a gain at the frequency w,, £ is a net-loss 
coefficient and a œ g (wp) /5Q? takes into account a frequency-dependence of gain coefficient in 
the vicinity of wy defined by a gain bandwidth ôQ. In such an oversimplified model, only one 
mode with the maximum net-gain o =g(w > w,)-€=0 is generated because the gain coefficient 
is energy-dependent, that is, it decreases with A (i.e. a gain is saturable that results in a mode 
competition or mode selection, Figure 1). 


Figure 1. Comb of frequencies generated by a laser with the repetition rate ôw. The gain band is centred at wọ. Mode- 
locking consists in excitation and synchronisation of w + ndw sidebands. 


In the case of active mode-locking, a periodic external modulation with the frequency of dw 
excites the +dw sidebands for each mode in the comb so that the modes A(w), A(w +6w) become 
coupled. In the framework of our oversimplified model, a steady-state regime AA = 0 is 
described by the equation in time domain [59-61]: 


(8-0) A(t) +a =z- vt A(t) = 0, (1) 


which is the classical equation for an oscillator in the potential defined by v œ 6w*. This 
equation has a trivial solution in the form of a Gaussian pulse [59—63]: 
A(t)= A expl -21n(2) t?/ 77), where the pulse amplitude A, is defined by the condition of energy 
balance of ø œ vt? (the saturable gain coefficient is energy, i.e. Af -dependent) and the pulse 
width t ova /v. In the general case, the excitation of A(t) in the form of Hermitian—Gaussian 
solutions of Equation (1) is possible. Since the pulse width is defined by v so that 
t <1/6Q6w, the minimum pulse widths of ~1/5Q are hardly reachable due to limitation 
imposed by 6w-value. The situation can be changed in the presence of the self-phase modula- 
tion (SPM) [65] and the dynamic gain saturation. Then t œ 1 [50 [66]. 


167 


168 


Fiber Laser 


Using the non-linear processes such as SPM, loss and gain saturation allows generating the 
ultrashort pulses due to mechanism of the so-called passive mode-locking [60]. Periodical 
perturbations caused by transitions through non-linear laser elements such as saturable 
absorber or gain medium enrich the spectrum with new components wọ + môw (m = 1, 2,..., M), 


which becomes locked through non-linear interaction [64]: 


0A(Q, ,Z i © œ œ% 
L a Zen(o feat dodo 1 X’ (0,0,0,0) 


g (2) 
xA(@,,z) A(@,, Ne ze "eelo -o - @,- 0), 


where a four-wave non-linear process defined by non-linear susceptibility x? mixes the 
frequencies w,, w, w; and w, during the propagation through a non-linear medium along the 
z-coordinate (Ak =k(w,)—k(w,)-k(w;)-k(w,) is the difference of wavenumbers, c is the speed of 


light, n is the frequency-dependent refractive index). 


Both active and passive mode-locking concepts can be easily united from the point of view of 
spice-time duality [64,67—69]. For instance, let’s consider heat diffusion equation: 


2 
See ous Soo -vxu, (3) 


where a heat radiated at x = 0 by the point source o diffuses along x-axis and is absorbed by 
cooler with the parabolic ‘cooling potential’. The replacements t > z and x > t result in an 
equation for ‘diffusion’ of light describing an active amplitude mode-locking (see Eq. (1)): 


A PA, (4) 


JÄ | EA A X 
Eq. (4) is clearly understandable in the Fourier domain: 57 =0A+v-—z-aw ?A, where an 
Z Ow 


external (i.e. active) modulation defined by the v-coefficient ‘diffuses’ (i.e. broadens) a field 
spectrum A and such diffusion is compensated by spectral dissipation defined by the a- 
coefficient. In the time-domain, the spectrum broadening corresponds to a light pulse short- 
ening due to parabolic potential action which is balanced by spectral dissipation causing a 
pulse widening. 


The space-time duality can be extended further with the help of diffraction-dispersion duality: 


ðA i @A ðA i 
= EA 
Oz 2k ax? Oz 


a / (5) 
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where k and £, are the wave number and group-delay dispersion coefficients, respectively. 
Both processes describe the beam/pulse spreading with propagation which is accompanied by 


phase ¢ profile distortion, i.e. by appearance of the chirp Q «dp /dx? (or d*/dt*). The active 
phase modulation from this point of view 


2 a -ivť?A (6) 


> a= 
(compare with (4)) looking as SL oA tiv 5-00 *A in the Fourier domain describes a 
‘diffraction’ (dispersion) in the frequency domain inspired by phase modulator which is 
balanced by spectral dissipation. The main difference from (4) is that the phase modulation in 
(6) distorts the phase and thereby produces chirp like the action of thin lens but in the time 
domain. In other words, the phase modulation in (6) pushes the spectral components out of 
the point of stationary phase 0¢ /dt=0, adding the frequency shift œ vt (Doppler shift) which 
enhances the spectral dissipation on the pulse wings and thereby forms a pulse like the active 
amplitude modulator. But the phase profile @(t) is parabolic in this case. 


The transition to a passive mode-locking looks straightforward, but one has to be careful in this 
case. The spice-time duality suggests a simple way to realize the temporal focusing like that 
in space domain: combination of phase modulation (‘time lens’) from Eq. (6) with dispersion 
(‘time diffraction’) from Eq. (5) allows compressing a pulse. Therefore, a replacement of time 
focusing by a time self-focusing (SPM) would provide a laser pulse self-trapping like the effect 
of laser beam self-trapping: 


ôA i, A 


= r ar iy|Al A, (7) 


which is the famous non-linear Schrodinger equation describing propagation of optical solitons 
in a fiber (6, < 0 corresponds to an anomalous dispersion, y is a SPM-coefficient) [35,70,71]. 


It is appropriate to mention here that the space-time duality x > t allows extending the 
physical context of consideration beyond scopes of optics. For instance, A, E = fdx |A|*, and 
p can be related for a mean-field amplitude, number of particle (mass of condensate) and 
momentum (wave number) for a Bose-Einstein condensate [39]. Then, it is clear that the 
dispersion and SPM terms in Eq. (7) describe the kinetic energy and four-particle interaction 
potential for gas of bosons. Such an interpretation opens a road to a quantum theory of solitons 
[72-74]. 


Following the same procedure for Eq. (4), describing the active amplitude mode-locking results 
in the simplest version of equation for a passive mode-locking, so-called cubic non-linear 
Ginzburg-Landau equation [35,36,75]: 
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2 
E ear A 


> = +K|Al A. (8) 


This equation describes a combined action of saturated net-gain (o), spectral dissipation (a) 
and non-linear gain (x). The last term results from loss saturation in a non-linear absorber with 
the response time much lesser than the pulse width. As will be shown below, such an as- 
sumption is valid for a broad class of fibre mode-locking mechanisms. Physics of passive mode- 
locking resembles that of active one: self-focusing in time domain causes a spectrum 
broadening which is balanced by spectral dissipation. Loss and energy-dependent gain are 
required for developing and stabilizing the mode-locking (all these factors are included in o- 
term which is < 0 for a steady-state pulse). Eqs. (7) and (8) have a similar solution 
A(t) oc sech(t /t) but the mathematical structures of these equations differ substantially that 
created discrepancies between concepts of the ‘true’ [76] and dissipative solitons (DSs, see next 
section) [36]. Combining Eqs. (7) and (8) gives the famous complex cubic non-linear Ginzburg- 
Landau equation (cubic CNGLE) [35-37,42]: 


OA 
+x Al A, (9) 


ôA -(4 OA 


gees jaf a jroasa 


which isa playground for study of DSs. Equation (9) allows a number of further generalizations 
such as: (i) description of non-distributed evolution due to dependence of the equation 
coefficients on z [77]; (ii) generalization of non-linearity type aimed first of all to adequate 
description of different mode-locking mechanisms (see below); (iii) taking into account the 
higher-order dispersions, i.e. w-dependence of f, [68]; (iv) taking into account the vector nature 
of light, i.e. transition to a system of coupled two-component CNGLEs [68,78-81], etc. 


Now let us consider the mode-locking mechanisms for fibre lasers in more detail. Active 
mode-locking can be utilized for DS generation from a fibre laser [82-84], but the wide- 
spread mechanism is based on the non-linear polarization rotation (NPR) which uses the effect 
of intensity-dependent polarization mode coupling in a fibre [85-88]. There is volumi- 
nous literature concerning the experimental realization of NPR mode-locking in fibre lasers; 
therefore, our selection of references is rather subjective and concerns the DS context [95- 
111]. 


It is known [70] that an ideal single-mode fibre supports two degenerate orthogonally 
polarized modes. However, a real fibre has inherent birefringence caused by core asymmetry 
or mechanical stress (Figure 2). 


Since SPM as well as cross-phase modulation (XPM) contribute to refractivity index with the 
strength defined by field intensity, such a contribution will change the state of polarization 
(SOP, Figure 3) [60,70,89] that can be described by coupled equations for two orthogonal (x 
and y) polarization components [70]: 
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Noncircular core Mechanical stress 


Figure 2. Fibre cross-sections with cylindrical symmetry broken by manufacturing process that leads to fibre birefrin- 
gence. 


polarizer fiber <0; analyzer 


Dw 9B 


wave plate SOP 


Figure 3. Block scheme of a self-amplitude modulator (SAM) based on NPR. Both linear birefringence and NPR con- 
tribute to the change of state of polarization (SOP, rotation angle 0) that leads to intensity-dependent transmission of 
this scheme. 
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where the dissipative factors from Eq. (9) are taken into account and AB=2r /L , describes a 
‘strength’ of linear birefringence (L, is a beat-length). As was shown in [81,91-94,179], the multi- 
scale averaging technique allows reducing Eq. (10) to the modified scalar non-linear Ginzburg- 
Landau equation (so-called sinusoidal Ginzburg-Landau equation) in which the self-amplitude 
modulation term (SAM, last term in Eq. (9)) is replaced by oc logQ( | A|?)A, where Q is a 


complex function defined by birefringence and settings of laser wave plates and polarizer. 
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Such an approach opens a way to multi-parametrical optimization of fibre lasers mode-locked 
by NPR. 


Despite its relative simplicity in principle as well as possibility of all-fibre-integrity of a laser, 
NPR in the form presented in Figure 3 is too sensitive to laser setup, uncontrollable perturba- 
tions and requires a precise manual tuning. The modified SAM setup, which can utilize both 
NPR and scalar SPM, is shown in Figure 4. It is the so-called non-linear optical loop mirror or 
figure eight laser (Figure 4) [112,113,265]. In principle, this setup is an all-fibre realization of 
additive-mode-locking [82,114] with inherently adjusted linear optical propagation lengths for 
counter-propagating beams. The main control parameter here is the beam splitting ratio p 
controlling the mutual intensities of counter-propagating beams. 


The unique property of this SAM setup is its ability to utilize different types of non-linearities 
for mode-locking (e.g. see [115-118]). Different modifications of this mode-locking mechanism 
have been used in DS fibre lasers [119-125]. Nevertheless, a fibre loop defining SAM remains 
environment- and tuning-sensitive. 


laser input/output 


Figure 4. Block-scheme of SAM based on a non-linear optical loop mirror. Splitter splits input laser beam into two 
counter-propagating ones (red and green arrows) with some splitting ratio p. The beams interfere after round-trip and 
partially return into a laser. The result of interference is intensity-dependent due to NPR or/and SPM within a loop. 


There is a class of alternative approaches utilizing non-fibre well-controllable non-linearities 
for mode-locking by the cost of broken fibre-integrity of a laser. Such an alternative was 
provided by development of high-non-linear semiconductor saturable absorber mirror (SESAM) 
[126-135]. The point is to put a semiconductor layer into a composed multi-layer mirror with 
the well-controllable spectral characteristics as well as with the adjustable intensity concerta- 
tion of penetrating field within a semiconductor layer. In fact, it is an advanced non-linear 
Fabri—Perot interferometer with the reflectivity coefficient depending on the incident intensity 
(or energy) [136]. Interaction of light with a semiconductor layer can be characterized roughly 
as excitation of carriers from a valence band of semiconductor to its conduction band. Excited 
carriers thermalize inside a conduction band with the character time ~100 femtoseconds. This 
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time defines a fastest response of SESAM to a laser radiation. Then, the thermalized carriers 
can relax into valence band or intra-band trapping states with the characteristic times from 
picoseconds to nanoseconds. Thus, SAM due to SESAM is slow in comparison with that due 
to NPR or SPM because the response times of the lasts are defined by intra-atom polarization 
dynamics, i.e. these times belong to femtosecond diapason. Additionally, the spectral diapa- 
son of SESAM response is substantially squeezed in comparison with that of pure electronic 
non-linearities due to resonant character of SESAM non-linearity. This can trouble the mode- 
locking within a spectral range exceeding the SESAM bandwidth. But the reverse side of the 
SESAM-band squeezing is that a non-linear response of SESAM becomes resonantly enhanced. 
This means that SESAM can provide more easily starting, stable and controllable mode- 
locking. The key characteristics of SESAM are [127] loss saturation fluence E,=hv /20, (hv isa 
photon energy, o is an absorption cross section), modulation depth ug=0,N (N is a density of 
states in semiconductor), relaxation (recovery) time T, unsaturable loss, saturable loss bandwidth 
and level of two-photon absorption. 


Akin mode-locking methods providing full fibre-integrity, broadband absorption, sub- 
picosecond response time and avoiding a complex multi-layer mirror weaving use nanotube 
and graphene saturable absorbers [30,137-143] and other low-dimensional structures [144]. 


From the theoretical point of view, the response of saturable absorber (SESAM or other 
quantum-size structures) to a laser field can be very complicate. In principle, one has to take 
into account finite loss bandwidth, its dispersion, dependence of refractive index on carrier’s 
(or exciton’s) density (so-called linewidth enhancement), complex kinetics of excitation and 
relaxation, etc. However, the praxis demonstrated that a simple model of two-level absorber 
is well working [145]: 


= A 
du UM fy ul (11) 


with some possible modifications (e.g. see [146]). Since DSs, as a rule, have over-picosecond 
widths (see next section), one may use an adiabatic approximation for (11) so that the expres- 
sion for SAM coefficient in the last term in Eq. (9) has to be replaced: 


Ho 


=) So 0 = 4 
il a 1+|al’ /P. 


(12) 


where P,=E,/T, is a loss saturation power. 


One may propose a hypothesis that an analogue of Kerr-lens mode-locking, which is a basic 
mechanism for generation of femtosecond pulses from solid-state lasers [60,85,147], can be 
realized in a fibre laser as well. Such an insight is based on possible enhancement of the laser 
beam spatial-trapping induced by non-linearity in a medium with spatially inhomogeneous 
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gain/loss or refractivity [148-152]. The model for analysis of such phenomena can be based on 
extension of dimensionality of Eq. (9), with taking into account the diffraction and transverse 
inhomogeneity of gain, loss or/and refractive index (the last can work as SAM due to the 
waveguide leaking loss) [153]: 


ôA -(4 CA 108A 


APA A CA 2 
az 262 kee a Ale BETSAN a 


t 


where cylindrical symmetry is assumed, x is a radial coordinate and x is a coefficient (complex 
in general case) which describes a transverse inhomogeneity of a fibre. Figure 5 shows the net- 
gain profiles (a) and the intra-laser pulse energies (b) as function of an effective aperture size 
obtained on the basis of variational approach for Eq. (13) [153]. The results demonstrate a 
principal feasibility of the Kerr-lens mode-locking regime for a DS fibre laser. 
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Figure 5. (a) Transverse net-gain profiles for different transverse parabolic distributions of net-gain coefficient which 
can be realized by inhomogeneous doping of fiber or by impact of waveguide leaking loss. (b) The dependence of in- 
tra-laser DS energy for a DS Yb-fibre laser with the 14 nm filter bandwidth and the average GDD of 330 fs?/cm [153]. 
DS collapses for large energies and cannot start for large aperture sizes (here SAM has inverse sign and the continu- 
ous-wave generation prevails). 


All these mode-locking techniques are realizable for both soliton proper and DS fibre lasers 
(excluding the Kerr-lens mode-locking which requires sufficiently high pulse energies 
provided by only a DS laser). Now let's consider the DSs fibre lasers proper. 


3. DS concept: Theory and experiment 


A ‘classical’ soliton can be formally defined as a solution of non-linear evolution equation 
belonging to discrete spectrum of the inverse scattering transform [71,76,154]. The non-linear 
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equations, which can be solved by inverse scattering transform, are ‘exactly integrable’. This 
means that they are akin to linear equations in some sense. In particular, they obey the 
superposition principle and, as a result, can be canonically quantized [155,156]. One has to 
note that integrability of a non-linear evolution equation and non-dissipative (non-Hamilto- 
nian) character of the latter are not equivalent because there are both non-integrable Hamil- 
tonian systems and integrable dissipative ones [36]. The point is that the DS concept is not 
connected with ‘integrability’; therefore, DSs are not ‘true’ solitons in a mathematical sense. 
However, many properties of DSs, in particular, their stable localization, robustness in the 
processes of scattering and interaction, well-organized internal structure, etc., resemble the 
properties of ‘true’ solitons. Formally, one may define DS as a localized and stable structure 
emergent in a non-linear dissipative system far from the thermodynamic equilibrium [36]. DSs are 
abundant in the different natural systems ranging from optics and condensed-matter physics 
to biology and medicine. In this sense, one may paraphrase that DSs “are around us. In the 
true sense of the word they are absolutely everywhere” [157]. Therefore, the concept of DS 
became well established in the last decade [36,37,42,158]. 


Stability of a DS under condition of strong non-equilibrium can be achieved only due to well- 
organized energy exchange with environment and subsequent energy redistribution within a 
DS. It results in energy flux inside a DS and, thereby, in DS phase inhomogeneity [36]. For a 
simplest case of Eq. (9), which has a DS solution in the form of A(t)=Ajsech(t / t)" 


(Y < Qt2=170?/dt7 is a dimensionless chirp parameter) [85,159,160], the DS energy genera- 
tion [36] 


oA 
ot 


2 2] ql? 
E=2o0|Al +2«|A]" 2a ra (14) 


as well as the spectrum | A(w) |? are shown in Figure 6 in dependence of f, (the data are based 
on an approach of [160]). One can see that the spectrum broadening transfers the action of 
spectral dissipation on a pulse ‘in whole’ into well-structured energy exchange: inflow at pulse 
centrum and outflow on its wings. The key characteristic of a dissipation inhomogeneity is a 
chirp, i.e. an inhomogeneity of phase. In absence of the chirp, the spectral dissipation acts on 
a pulse in whole that, in particular, induces a multi-pulse instability [161]. However, a power- 
dependent chirp causes inhomogeneity of energy transfer (Figure 7). Energy flows in the 
region closer to central wavelength where the gain is maximal. This region is located in the 
vicinity of pulse maximum. Energy flows out from the spectrum wings which are located on 
the wings of pulse, that is, the pulse localization is supported by spectral dissipation through 
non-linear mechanism of chirping [42,160,162]. One has to note that a direction of energy fluxes 
inside a DS depends on parameters and can be inversely related to the direction shown in 
Figures 6,7 (i.e. energy can flow from wings to centre). The corresponding structure was named 


dissipative anti-soliton [210]. 
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Thus, an additional mechanism of SAM (in addition to mechanisms considered in the previous 


section) appears, which provides unique robustness of DSs (i.e. DS exists within a broad range 
of laser parameters [163,164]). 


(b) 


arb. un. 


Log(A()). 


Figure 6. (a) Profile of energy generation and (b) logarithm of spectral power in dependence on GDD (normal disper 
sion range) for a DS of [160]. 


Below, we will consider a chirp as the essential characteristic of DS [210]. One of the reasons 
is that the chirp allows DS to accumulate energy E œ, which means that DS is energy- 
scalable [37,164-167]. The last statement does not mean that a chirp-free pulse is not energy- 
scalable. However, the energy-scalability of such pulses can be provided by only fine- 
tuned and separated control of SPM and GDD that can be achieved in solid-state oscillators 
[26,168] or in large mode area (LMA) fibre lasers [169]. For fibre lasers such an approach 
entails issues of full-fibre integrability, higher-order mode control [170,171]* and thermo- 
effects impact [172]. 


© 


A(@); arb. un, 


Figure 7. (a) Profile of energy generation and (b) power in dependence on chirp for a DS of Eq. (9): o = —0.01, x = 0.025, 
a=0.05, A, =1,t=1. 


2 However, namely LMA and photonic-crystal fibres could realize a Kerr-lens mode-locking in a fibre laser [152,153]. 
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In the terms of space-time duality (see above), the mechanism of formation of time window, 
within which a DS is localized, resembles a phenomenon of total internal reflection from some 
‘borders’ created by phase discontinuity. Such borders are formally defined by the equivalence 
of the wave number of out-/in-going radiation k(w)=B,w 7/2 (wave number of dispersive linear 
wave) and the DS wave number q=y P) : k(+A)=q, where P)= | Anax |? is a DS peak power and 
a DS spectral width is A=,/2yP)/f, [173]. Since a system is dissipative, the above phase 


equilibrium has to be supplemented by loss compensation condition: spectral loss œ aA? has 
to be compensated by non-linear gain xP). Combination of above criteria gives a definition of 
the parametric limits for DS [44,173]: 


2a 2/3, E> 
aar cfa , (15) 


Bx | 2,E>0 


where E is a DS energy. Eq. (15) is valid for the cubic-quintic CNGLE in which SAM has a form 
of (x | Al?-x¢ | Al4)A, that is, a non-linear gain is saturable (despite the unsaturable SAM in 
Eq. (9)). A SAM saturability is necessary for DS stabilization [174]. Equality in Eq. (15) 
corresponds to the DS stability border where o = 0 (see Eq. (9)). The asymptotic E > 
corresponds to a perfectly energy-scalable DS or to a phenomenon of dissipative soliton reso- 
nance (DSR) [37,44,165-167,175-183], which is sufficiently robust, exists in different SAM 
environments and even within the anomalous GDD range [178,182]. Important property of 
DSR is that the DS energy E can be scaled without loss of stability by plain scaling of laser 
average power or/and its length L [44,180,184]. The chirp scales with length as well. As a result, 
the DS peak power and spectrum width tend to a constant for fixed parameters of Eq. (9) (i.e. 
fixed a, f, o and x) and the energy scaling is provided by DS stretching in time domain. 


This ideology of energy scaling by the pulse stretching goes back to the so-called wave- 
breaking-free or stretched pulse fibre lasers where the propagation within the anomalous- 
dispersion fibre sectors alternates with the propagation under normal GDD action [96,101,185— 
187]. As a result of pulse stretching, the non-linear effects in such systems are reduced, which 
allows increasing an energy and suppressing a noise. As an alternative approach, one can 
exclude an anomalous GDD at all and to realize a so-called similariton regime, when a pulse 
accumulates an extremely large chirp and, thereby, an energy [33,103,188]. However, a self- 
similar regime is not soliton-like one in nature; therefore, we will focus on the all-normal- 
dispersion fibre lasers (ANDi) which produces DSs possessing a high stability within a broad 
range of laser parameters [97,101,185,186,189,190]. In Figure 8, the energy-scalable DS lasers 
are sub-divided into three main types: (1) all-fibre, (2) fibre with a free-space sector and (2) LMA 
including rod and photonic-crystal fibre PCF. 


The advantage of the first type of lasers is their integrity, which does not require an operational 
alignment, includes potentially compressing and delivering sections, and provides environ- 
ment insensitivity and easy integrability with fibre-amplifier cascades [200]. The last advant- 
age is especially attractive because it allows a direct seeding of DS into chirp-pulse amplifier 
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without preliminary pulse stretching. The second type of the DS fibre lasers can be considered 
as a testbed for development of the first type. No wonder that the results achieved here are 
more impressive (Figure 8). At last, the third type of the DS fibre lasers is most akin to the thin- 
disk solid-state ones with simultaneous advantage of the broad gainbands. Such lasers provide 
the DS energy scalability by scaling of laser beam area in combination with the scaling of laser 
period and average power. Nevertheless, one has to keep in mind that both LMA and PCF 
technologies have some disadvantages (see above) which make them similar to solid-state 
lasers. 
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Figure 8. Experimental realizations of DSs in fibre lasers (Refs. 
[7-9,82,83,100,102,105-107,109-111,122,124,125,130,132,135,138,139,140,141,187,191-209,232]; Ref. [111] corresponds to a 
DS in the anomalous GDD region). 


The diversity of the results obtained (Figure 8) needs a comprehension from a unified view- 
point; therefore, let’s survey briefly some theoretical aspects relevant to the DS fibre lasers. 
There is vast literature regarding the theory of DSs. Some preliminary systematization can be 
found in [44]. However, it is necessary at first to declare the stumbling block of this theory: 
absence of a unified viewpoint. There exist unbroken walls between the circles of scientific 
community exploiting and exploring the DS concept: walls between the solid-state and fibre 
laser representations of the theory, condensed-matter one, numerical and analytical ap- 
proaches, etc. Briefly and conditionally, the relevant theoretical approaches can be divided 
into (1) numerical, (2) exact analytical and (3) approximated analytical. The last includes the 
models based on (3.1) perturbative, (3.2) adiabatic models (AM) as well as those based on (3.3) 
phase-space truncation (i.e. variational approximation (VA) and method of moments (MM)). 
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As was emphasized repeatedly, both linear and non-linear dissipations are crucial for the DS 
formation. The simplest and most studied models for such a type of phenomena are based on 
the different versions of CNGLE (e.g. Eq. (9)). 


Extensive numerical study of DSs of the cubic-quintic CNGLE has been carried out by N. N. 
Akhmediev with co-authors [35-37,42,157,166,175,176,178,182]. The simulations have allowed 
finding the DS stability regions for some two-dimensional projections of CNGLE parametrical 
space. The summarizing description of the results obtained is presented in [44]. Most impres- 
sive results are: (i) parametric space of DS has a reduced dimensionality resulting, in particular, 
in the appearance of DSR; (ii) DSR remains in a model with lumped evolution that is typical 
for the most of fibre lasers; (iii) DSR and, correspondingly, DS exist within the anomalous GDD 
region as well. However, the main shortcomings of the numerical approaches are: the para- 
metrical space under consideration is not physically relevant, and the true dimensionality of 
DS parametric space is not identified. It is clear that the only advanced and self-consistent 
analytical theory of DS would provide, in particular, a true representation of DS parametric 
space and DSR conditions. 


As was mentioned above, the evolution equations describing DSs are not-integrable. The 
efforts based on the algebraic techniques [62,213,214] and aimed to finding the generalized DS 
solutions of CNGLE were not successful to date. Nevertheless, few exact partial DS-solutions 
are known. For instance, sole known exact analytical DS-solution of cubic-quintic CNGLE is 
[110,166,176,182,189,211,212]: 


Ap ; 
m cosh(t/r) zB op In(cosh(t/r) + B) + a (16) 


where A, B, t, Y and ¢ are real constants [189]. This solution belongs to a fixed-point solution 
class, which means that it exists only if some constraints are imposed on the cubic-quintic 
CNGLE parameters. Solution (16) provides with important insights into properties of DSs. In 
particular, the systematical classification of DS spectra (truncated concave, convex, Lorentzian 
and structured spectra) and DS temporal profiles (from sech-shaped A(t) to tabletop one) is 
possible in the framework of analytical approach. The transition to a DSR-regime reveals itself 
in the ‘time-spectral’ duality shown in Figure 9 [141,189]. The sense of this ‘duality’ will be 
explained below from the point of view of adiabatic theory of DSs. 


The crucial shortcoming of the approach based on few exact DS solutions of evolution 
equations is that the strict restrictions are imposed on the equation parameters. As a result, the 
DS cannot be traced within a broad multidimensional parametric range and the picture 
obtained is rather sporadic and is of interest only in the close relation with the numerical results 
and experiment. Some additional information can be obtained on the basis of perturbation 
theory which provides with a quite accurate approximation for a low-energy DS [215-217]. 


Most powerful approaches to the theory of DSs have been developed in the framework of 
approximated techniques (for review see [44]): AM [165,167,217-221], VA [77,177,222-224,225] 
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Figure 9. ‘Time-spectral duality’ for an energy-scalable DS: bell-like time-profile (blue curve on the left picture) trans- 
forms into tabletop one with truncated edges (right) and, inversely, truncated tabletop spectrum (red curve on the left 
picture) transforms into bell-like one (right) with the DS energy growth. 


and MM [175,210,222,226]. The most impressive results obtained are: (i) physically relevant 
representation of DS parametric space was revealed (it is a so-called master diagram, see [44] 
for review and Figure 10); (ii) such a representation allows understanding the structural 
properties of DS and its energy-scaling laws (i.e. DSR conditions) for different mode-locking 
techniques; (iii) DS dynamics and an issue of optimal arrangement of laser elements providing 
the maximum DS stability and energy have been explored [224,225,227,228]; (iv) vectorial 
extension of VA concerning a vector DS (VDS) was endeavoured [229]. 
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Figure 10. DS master diagrams for the cubic-quintic CNGLE (black) and the CNGLE with a SAM defined by Eq. (12) 
(red) (in the last case K =1 I yP ;). DSR ranges correspond to a so-called positive branch of DS [44,167], which has a 
highest DS energy-scalability and stability [221]. The dashed curve corresponds to a DS stability border obtained from 
numerical simulations of cubic-quintic CNGLE taking into account a quantum noise [243]. Dot blue curve shows the 
DS border under effect of SRS [260]. The parametrical space shown is the physically relevant parametrical space of DS. 
Points correspond to different scenarios of DS destabilization (see text). 
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Both AM and VA demonstrate two-dimensional representation of DS parametric space in the 
form of master diagram. Dimensionality can grow with complication of CNGLE non-linearity 
when SPM becomes saturable so that the cubic non-linear term in Eq. (9) has to be replaced by 
[(x-x¢ | Al?)-ily +x I A1?)]I A l?A [217]. This effect can appear ina fibre laser with NPR (e.g. 
see [179] where such a completely cubic-quintic CNGLE is connected with the NPR mode- 
locking technique). In this case, DS soliton exists in both normal and anomalous GDD regions 
[175,182,233]. 


The master diagram is a manifold of isogains (i.e. curves with o =const <0). Figure 10 demon- 
strates the zero-isogains (ø = 0) corresponding to the DS stability limit (upper curves) as well 
as the borders between ‘energy-scalable’ and ‘energy-non-scalabe’ branches of DS (lower 
curves, see [44,167] for a formal definition’). Figure 10 demonstrates that the saturation of SAM 
(so-called reverse saturable absorber provided, for instance, by NPR or graphene [141]; black 
curves) enhances the DS stability in comparison with an unsaturable SAM (red curves). Since 

lim E=% (const € the DSR range and its maximum value are of 2/3) for a saturable SAM (cubic- 


C-const 

quintic CNGLE, black curves in Figure 10). Such a property of isogain curves corresponds to 
the DSR phenomenon. Since the stability threshold is defined by the condition of 
C=2ay | fx <2/3, one may conclude that the broadening of spectral filter band (or gainband) 
enhances stability against multi-pulsing (a 1/507, see above) [107,108]. Simultaneously, 
SPM has to be balanced by GDD (C œ y /£,) that, in combination with normalization of energy, 
gives the energy-scaling law E œ L along a DSR curve. VA predicts [230]: 


Ex B,/Va (17) 


which agrees with experimental observations of linear growth of DS energy with bandwidth 
[107,108] as well with a rule Eœ L since B,«L in the frameworks of distributed CNGLE. 


In the case of unsaturable SAM corresponding to SESAM, some nanotube and graphene 
absorbers, Kerr-lensing, etc. (see Eq. (12)), the energy scaling requires scaling of the control 


parameter C. In this case, the asymptotic energy scaling law for Ex 2 ya > 1 becomes [44]: 
Ex B,/CVa. (18) 


The spectral properties of DS are described clearly in the frameworks of AM [44,167]. In the 
simplest case of cubic-quintic CNGLE, the DS spectrum A(w) in the limit of |y | >1 is a 
Lorentzian profile which has a characteristic width Q, and is truncated at frequencies +A 
[44,167,218]: 


P H A? — 2 
oe (19) 


2 2 
Lt 


3 Energy-non-scalable branch has two distinguishing characteristics: it turns into solution of Eq. (9) with C, 7-0 
(‘Schrödinger limit’ [218]) and is unstable in absence of dynamic gain saturation, i.e. if o is not energy-dependent [221]. 
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where H is a Heaviside function. The DS energy is 
Eœ arctan (4/9; )/Q, ; (20) 


and 


2l dt 
-Anaa 8) ey 


Here, we trace the zero-isogain o = 0. The DS time-profile is defined by an implicit expression: 


t= arctan $2 + A arctan S84) F (22) 
A Q A 


L 
with the DS width of t œ AHA? +0? )!. Now, there are the following limiting cases: 


E>0 
P, >00 
_ (aed 
a Q, — const i (23) 
T>0 


y|AP =q- BA? tanh? (t/z) 


It is clear that in this ‘low-energy’ sector the DS time-profile is bell-like and its spectrum has 
tabletop form (Q, > A). In the DSR limit, one has: 


E > œ 

P, >1/¢ 

lim — 71 BC 24 
c32/0, 30 i (e 
T — const 


y|Al =q — pO? tan? (tQ, /14) 
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that is, a DS in the DSR sector has a flattop temporal profile and a Lorenzian spectrum 
(Q, <A). Eqs. (24) demonstrate that asymptotical growth of DS energy leads to a spectral 
condensation (Q, — 0) without a parallel temporal thermalization (t °°), which means an 
inevitable destabilization of a plain energy-scalability [180]. This conclusion does not mean a 
participial impossibility of DS energy-scaling in the frameworks of cubic-quintic CNGLE 
model. For instance, a saturable SPM allows DSs with tabletop profiles and | dọ /dt | > œ on 
the pulse edges. Such a DS possesses enhanced energy scalability and was observed experi- 
mentally [231]. 


4. DS spectrum and stability 


As was explained, the dual balances in frequency domain: 


ee = (25) 


ad? <> KP, 


are formative for DS existence and stabilization. No wonder that the spectrum of DS is 
benchmark of its inherent properties. 


Figure 11. Perturbed DS spectrum [241]. 


Prior to consider the aspects of interweaving of spectral and stability properties of DSs, one 
has to point to a possibility of multi-wavelength multi-pulsing DSs provided by DS robustness. 
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As was demonstrated in [234] theoretically, the multi-DSs compounds in a mode-locked laser 
can be stabilized at multiple frequencies. Experimentally, such multi-frequency DS com- 
pounds can be realized by birefringence filters with a periodical (interference-like) dependence 
of transmission on wavelength under conditions of sufficiently broad gainband and powerful 
pump [235-239]*. 


Figure 12. Exploding DS corresponding to parameters of point A in Figure 10. Left: contour-plot of instant power, 
right: 3D-graph of instant power in dependence on local time t and propagation distance z (arbitrary units) [243]. 


As was demonstrated in previous section, DS has non-trivial internal structure due to energy 
fluxes inside it. The elements of this structure (internal modes) can be excited that causes 
pulsating or chaotic dynamics of DS with preservation of its temporal and spectral localization 
[241]. The spectral envelope acquires a shape of ‘glass with boiling water’ (Figure 11). Ap- 
pearance of such perturbations is understandable in frameworks of the DS perturbation theory 
in spectral domain [242]. One has to note that such perturbations take a place inside the DS 
stability region where o < 0 (below the corresponding upper curves in Figure 10). Above the 
stability boarder (‘no DS’ region in Figure 10), there are three main destabilization scenarios 
[243]. For small energies and in the vicinity of stability border (point A in Figure 10), DS is 
exploding (Figure 12), which means its aperiodic disappearance with excitation of continuous 
waves and subsequent DS recreation [37,244—248]. With the energy growth (point B in Figure 
10), the rogue DSs develop (Figure 13) [37,249-253]. Such a regime can be interpreted as DS 
structural chaotization, that is, generation of multiple DSs with strong interactions causing 
extreme dynamics. 


For sufficiently large energies in the vicinity of stability border (point Cin Figure 10), the typical 
destabilization scenario is the generation of multiple DSs (Figure 14). The source of this 
destabilization is the growth of spectral dissipation caused by DS spectral broadening with 
approaching to stability border so that the DS splitting becomes more energy advantageous 
[161]. Moreover, the DS splitting can be enhanced by its phase inhomogeneity because the gain 
(energy in-flow) is maximum at the points of stationary phase dọ / dt=0 [254]. Such a splitting 


4 A multi-porting configuration of a DS laser supports even simultaneous generation of conventional and dissipative 
wavelength-separated solitons [240]. 
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Figure 13. Rogue DSs corresponding to parameters of point B in Figure 10. Left: contour-plot of instant power, right: 
3D-graph of instant power in dependence on local time t and propagation distance z (arbitrary units) [243]. 
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Figure 14. DS molecule corresponding to parameters of point C in Figure 10. Left: contour-plot of instant power, right: 
3D-graph of instant power in dependence on local time t and propagation distance z (arbitrary units) [243]. 


can result in an extreme dynamics like solitonic turbulence (Figure 13) [255] or regular multi- 
DS complexes (Figure 14, so-called DS molecules) which can have non-extreme internal 
dynamics (soliton gas) and interact with a background (soliton liquid) [256-259]. If the dynamic 
gain saturation (see [44]) contributes, the multi-DS complexes can evolve slowly into a set of 
equidistant DSs with repetition rate multiple of laser one (harmonic mode-locking) [158,159]. 


The numerical simulations of cubic-quintic CNGLE with taking into account a quantum noise 
validated the fact of inconsistency of spectral condensation and absence of temporal thermo- 
lization that breaks the DS energy scalability (see previous section) [243]. As a result, the DS 
stability region breaks abruptly with energy growth (dashed curve in Figure 10) and multitude 
of turbulent scenarios of DSs evolution develops (Figure 15) [34,243]. 


Serious limitations on power and energy scalability of DSs in fibre lasers arise from stimulated 
Raman scattering (SRS) [2,109]. The stability border of DS under action of SRS is shown in 
Figure 10 by dot blue curve (DS is stable on the left of this curve) [260]. As was found, SRS 
enhances the tendency to multi-pulsing with energy growth caused by enhancement of 
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Figure 15. DS turbulent regimes (contour-plot of instant power, arbitrary units) corresponding to the parameters of 
points D, E, F, G and H in Figure 10 [243]. 


spectral dissipation due to SRS [260]. Simultaneously, generation of anti-Stokes radiation 
causes chaotization of DS dynamics and irregular modulation of DS temporal and spectral 
profiles [261] (Figure 16). DS profile remains localized, but it is strongly cut by colliding dark 
and grey soliton-like structures [34]. 


As was shown, the DS dynamics can be regularized by formation of dissipative Raman soliton 
(DRS). DRS can exist in the form of DS which is Stokes-shifted due to self-Raman scattering 
(Figure 17, a) [260] or as bound DS-DRS complex (Figure 17, b) [262]. In the last case, stabiliza- 
tion is achieved by feedback, i.e. reinjection of Stokes signal through a delay line [263]. In the 
absence of a feedback, the Raman pulse is noisy [264] (see dash line spectrum in Figure 16, b). 
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Figure 16. Wigner function of turbulent DS in the presence of SRS. 
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Figure 17. (a) Wigner function of single DRS [260] and (b) spectrum of bound DS-DRS complex [262]. 


5. Vector DSs 


As was pointed above, SOP can play leading role in a fibre laser dynamics. In particular, it can 
contribute to mode-locking or/and spectral filtering. However, diapason of polarization 
phenomena in a DS fibre laser spreads essentially broader. As was found, intrinsic fibre 
birefringence (Figure 2) can lead to DS splitting into two independent SOPs [78]. This phe- 
nomenon is used to realize the NPR mode-locking mechanism where a DS SOP evolves (or 
remains locked) as a whole during propagation [265-269]. The polarization dynamics can be 
fast (<T) or slow (>>T) and vary from regular (with possible period multiplication or harmonic 
mode-locking) to chaotic one [270,271]. There are evidences of ultrafast SOP evolution when 
SOP changes across a DS profile [272]. 


The specific multiple pulse instability of vector dissipative solitons (VDSs) leads to generation 
of the bound states of DSs with different SOPs (vector soliton molecules) which are locked by a 
non-linear coupling [273,274] or by a group-velocity locking produced by spectrum shift 
between DSs with different SOP [275]. As was shown experimentally (Figure 18) [276], the 
dynamics of VDS molecules can be highly non-trivial and demonstrate both fast and slow 
periodic switching between fixed SOPs as well as SOP procession, which is especially inter- 
esting for fibre laser telecommunications based on polarization multiplexing. 


The important breakthrough in the recent theory of VDSs is the demonstration of insufficiency 
of approaches based on the coupled CNGLEs (like (10)) for adequate description of DS 
polarization dynamics. It was demonstrated that an active medium polarizability contributes 
to DS dynamics substantially [277]. As was shown, the SOP-sensitive interaction between DS 
and a slowly relaxing active medium with taking into account the birefringence of fibre laser 
elements and light-induced anisotropy caused by elliptically polarised pump field change the 
SOP at a long time scale that results in fast and slowly evolving SOPs of VDSs (Figure 19). 


The non-trivial contribution of active medium kinetics and polarizability with taking into 
account the pump SOP and SPM demonstrates a complex dynamics including spiral attractors 
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Figure 18. Polarization dynamics of VDS molecules in an Er-fibre laser mode-locked by carbon nanotubes [276]. Top 
row demonstrates the evolution of Stokes parameters. The bottom row reproduces this evolution on the Poincaré 
sphere (each point corresponds to SOP after one laser round-trip). 


and dynamic chaos (Figure 20) [278]. One may assume that such a non-trivial polarization 
dynamics is of great importance for DS energy scaling, in particular, due to vector nature of 
SRS [279]. These topics remain unexplored to date. 
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Figure 19. Polarization dynamics of VDS with taking into account an Er-fibre polarizability [277]. Both fast and slow 
SOP dynamics exist in dependence on fiber laser birefringence strength. 
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Figure 20. Dynamics of active fibre inversion components, which are polarization-sensitive (n4 and n,,) and non-sensi- 
tive (no), with the corresponding slow SOP evolution of VDS [278]. 


6. Conclusion 


The recent progress in development of ultrafast fibre lasers and advances in exploring of DS 
are interrelated. DSs allowed scoring a great success in ultrashort pulse energy scalability that 
is defined by unprecedented stability and robustness of DS. At this moment, it is possible to 
achieve over-MW peak powers for sub-100 fs pulses directly from a fibre laser at over-MHz 
repetition rates. New spectral diapasons became reachable owing to development of mid-IR 
active fibres and using the frequency-conversion directly in a laser. Development of new mode- 
locking techniques, especially based on using of SESAMs, graphene and another quantum- 
sized structure allowed improving a laser stability, integrity and environment insensitivity. A 
great advance has been achieved in the theory of DSs. New powerful analytical techniques 
based on extensive numerical simulations and experimental advances extended understand- 
ing of the DS fundamental properties and revealed new prospects in improvement of charac- 
teristics of ultrafast fibre lasers. Based on achieved results, one may outline some unresolved 
problems. As was found, there are stability limits for a DS energy scaling imposed by optical 
turbulence and SRS. Deeper insight into the nature of these phenomena could allow to 
overcome these limits without substantial complication of laser setup. Simultaneously, control 
of intra-laser spectral conversion is a direct way to broadening of spectral range. Then, the 
dynamics and properties of VDSs remain scantily explored. Recent studies demonstrated a 
multitude of polarization phenomena, which cannot be grasped in frameworks of existing 
models. In particular, polarizability and kinetics of an active fibre in combination with 
birefringence of a laser in a whole can contribute non-trivially to a laser dynamics. As an 
additional aspect of further development, one may point at the development of new mode- 
locking techniques, which could improve DS stability and integrity of a fibre laser, decrease 
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pulse width and extend a diapason of pulse repetition rates. At last, one has to remember that 
a fibre laser is an ideal playground for study of complex non-linear phenomena and, undoubt- 
edly, new bridges between different fields of science will be built with a further progress of 


ultrafast fibre lasers. 
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